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Abstract: We study the distribution 1/a da/dQt, where Qt is the modulus of the trans- 
verse momentum vector, obtained by summing over all hadrons, in the current hemisphere 
of the DIS Breit frame. We resum the large logarithms in the small Qt region, to next-to- 
leading logarithmic accuracy, including the non-global logarithms involved. We point out 
that this observable is simply related to the Drell-Yan vector boson and predicted Higgs Qt 
spectra at hadron colliders. Comparing our predictions to existing HERA data thus ought 
to be a valuable source of information on the role or absence of small- a; (BFKL) effects, 
neglected in conventional resummations of such quantities. 
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1. Introduction 



With the imminent advent of the LHC, considerable effort is being dedicated to utiUsing 
existing colhder data and theoretical predictions for QCD observables to foster an even 
better knowledge of crucial parameters such as the strong coupling constant and parton 
distributions in hadrons Q. 

An invaluable source for such data is the HERA collider, which continues to play an 
important role in phenomenological studies of QCD. The recently concluded HERA-LHC 
workshop^ was in fact dedicated towards the aim of directly linking HERA QCD studies 
to those that will be important in a discovery context at the LHC. 

In the present paper we shall highlight one such study. To illustrate our point we choose 
the transverse momentum (Qt) spectrum of the Higgs boson for which very accurate theo- 
retical predictions exist in the literature |Q] . These estimates resum large logarithms in the 
small Qt limit to next-to-next-to leading logarithmic (NNLL) accuracy'^. The resummed 
predictions are then combined with fixed next-to-leading order (NLO) computations Q so 
as to achieve the best possible accuracy over a wide range of Qt ■ Such accurate studies are 
important in the context of formulating improved strategies to extract the signal and to 
enhance its statistical significance over backgrounds. 

One concern that has been a point for some discussion (see e.g. Ref. [Q]) is the role that 
might be played by neglected small-x effects that one may expect could be relevant at the x 
values involved in vector boson and Higgs production at the LHC. Here x can be interpreted, 
as usual, as the momentum fraction of the incoming hadron carried by the struck parton. If 
neglected small-x effects are indeed important in this context, they could in principle lead 
to broader Qt spectra than those predicted by conventional resummations (sometimes also 
referred to as Collins-Soper-Sterman or CSS formalism [^) alone. Henceforth we shall refer 
to these resummations involving Qt spectra in hadron-hadron collisions with the generic 
label Drell-Yan Qt resummations. 

It was observed for instance in Ref. [Q] that effects due to small-x enhancements that 
were suggested by phenomenological studies for semi-inclusive DIS (SIDIS) Qt distribu- 
tions, could be very significant (especially in the context of massive vector-boson produc- 
tion) when extrapolated to the smaller x values that will be important at the LHC. It was 
also suggested in Ref. |Q] that such effects could be visible with Tevatron (run 2) data if 
one concentrates on forward production of vector bosons alone rather than integrating over 
all rapidities. On the other hand, studies for many event-shape variables in the current 
hemisphere of the Breit frame in DIS have been very successful, down to moderately small 
X values {x ~ 10~^), based on resummation [0, ^, ^ that did not account for BFKL or 
formal small-x effects. The comparison to HERA data for several event-shape variables 
can be found in 

It is also clear however that event-shapes are somewhat different from Qt spectra of 
the Drell-Yan type, significantly due to their direct sensitivity to gluon emission, rather 

^bttp:/ /www. desy.de/~heralhc 

^For reference, most resummations e.g. those for several event-shapes at LEP and HERA, achieve only 
a next-to-leading (NLL) logarithmic accuracy. 
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than purely through recoil. This difference means that event-shape variables receive non- 
perturbative corrections that scale as 1/Q [|^, where Q is the DIS hard scale. These 
power corrections are due to soft gluon emission alone and hence are independent of x. 
The quantity we shall study in the present article is closer in nature to the Drell-Yan Qt 
spectra since the leading non-perturbative effects here scale as and can be associated 
to what is commonly known as "intrinsic /c^" of partons inside hadrons. Thus one would 
expect any missing small- a; effects that appear in the present context (perhaps as suggested 
in Ref . in terms of a small-x enhanced smearing of the conventional resummation) , to 
manifest themselves in a very similar way in the Drell-Yan case. It is thus conceivable that 
for the Qt distribution we present here, deviations are seen from the resummed form, even 
at X values already studied successfully for event-shape variables. 

To be more precise, the observable we study here is the distribution 1/ a da / dQt, with 
Qt being the modulus of the transverse momentum vector of all particles in the current 
hemisphere {T~ic) of the DIS Breit frame: 

4=J^4i> (1-1) 

where kt is measured with respect to the photon axis. Note that the addition over particles 
mentioned above, is vectorial in nature, in contrast to say that involved in event-shapes 
like the jet broadening observable 0, where one adds the modulii of individual particle 
transverse momenta. This is also a different quantity from the resummed z flow studied in 
Refs. [Q, ^ while also being directly related to the Drell-Yan Qt spectra and thus should 
provide an independent probe of small-x broadening effects. 

To relate the quantity here to Qt spectra in Drell-Yan type processes one notes that 
in the Drell-Yan case one has the massive vector-boson recoiling against emissions from 
two incoming hard partons, while for our DIS example, the transverse momentum of the 
current hemisphere is equal and opposite to that of the remnant hemisphere. If one assumes 
the remnant hemisphere Qt spectrum to be entirely generated by emissions from a single 
incoming leg the relationship to Drell-Yan Qt is obvious: we just have to account for the 
form-factor of one incoming parton instead of two and thus we have, at the level of form- 
factors, A^jg((5,Qt) = \J Apy(Q) Qt)- The variable indicates moment space, conjugate 
to X. 

The simple relationship we mentioned immediately above breaks down at the next-to- 
leading logarithmic (NLL) level, due to the non-global nature ^ of the DIS observable 
we have introduced. Non-globalness in the present case is a consequence of looking at just a 
single hemisphere and provides an extra factor in the DIS case, to do with soft emissions at 
large angles, which is absent for the global Drell-Yan quantity. We shall of course account 
for this factor, but we stress here that the fundamental relationship of our quantity to Drell- 
Yan Qt (which is exactly of the square-root form we wrote above, at leading-logarithmic 
accuracy) is unchanged by this complication. In particular any neglect of terms that are 
enhanced at small x, ought to be of similar significance in the two cases. 

To further this investigation we obtain here a resummed result for the above quantity, 
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to next-to-leading logarithmic (NLL) accuracy ^ and combine it with fixed-order predic- 
tions to 0(0^) accuracy. Our result is thus suitable for comparison with data over the entire 
measured range of Qt values. We also comment on the effect of smearing our result with 
a Gaussian function, as is usual to accommodate the so-called intrinsic kt of the incoming 
parton, which has a non-perturbative origin. The resultant prediction can then be directly 
compared to data which should be of interest especially at lower x values. If discrepancies 
are visible at low x then one may consider a small-x enhanced smearing function as was 
the case in Refs. [Q, This is certainly not a substitute for a detailed treatment based on 
a physical understanding of the small- j; region but merely a phenomenological investiga- 
tion into how such effects may be parameterised if present in the first place. Subsequently 
one may also consider the extrapolation of our findings to hadron colliders. We note that 
preliminary data from HI are already available and await the final versions together 
with potential data from the ZEUS collaboration. 

The outline of this paper is as follows: in section ^ we put together the different 
ingredients required to obtain our NLL resummed result which we compute in impact 
parameter or h space as is most convenient for resummations. Once we obtain the h 



space answer we find, in the following section, its transform to Qt space and comment on 



its main features. In section ^ we carry out the matching of our resummed result to the 



full O {of) result from the fixed-order program DISENT [14|. Finally we comment on the 
potential role of non-perturbative effects that are expected to take the form of a smearing 
of the Qt distribution with a function representing the "intrinsic transverse momentum" of 
partons inside hadrons. Here one may try different choices for the smearing function and 
search for any discrepancies at lower x values, between our results and the data. We shall 



leave the details of this to our forthcoming phenomenological investigation |15| . 
2. Resummation 

At Born level, the struck quark is aligned along the photon axis, the quantity in question 
{Qt) vanishes and the distribution is essentially a delta function: da/dQt ~ S{Qt). 

At small Qt the emission of soft and collinear gluons deform the delta function. One 
may, on general grounds, expect this deformation to take the form of a Sudakov form-factor. 
This is essentially true over a large range of Qt values with the caveat that at very small Qt, 
the correct result is no longer of Sudakov form. The reason for this is the Parisi-Petronzio 
observation that the smallest Qt values are in fact obtained by vectorial cancellation of 
emissions rather than by suppressing the transverse momenta of each individual emission 
p^. We shall explain this issue more quantitatively, with reference to our observable, in a 



later section . For now we proceed with resumming the large logarithms that arise at small 

Qt- 

To carry out the resummation, we have to address two distinct kinematic regimes: 

• collinear emissions (soft or hard) along the directions of the incoming and scattered 
quark directions. 



^Equivalently we seek single-loganthmic accuracy in the resummed exponent that we shall compute 
subsequently. 
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• soft emissions at large angles to the incoming and scattered (outgoing) quark axis 
(since in the Breit frame the incoming and scattered quarks are back-to-back at Born 
level) . This contribution arises due to the non-global nature [p!l| , |l^] of the observable 
and is a correlated multi-gluon emission term, that can only be computed in the large 
Nc limit. 

We shall treat each region in turn starting with the collinear enhanced contribution 
and then including the non-global term that arises from the piece of the fixed-order matrix 
elements with only soft enhancement (i.e. is integrable over soft gluon directions). 

2.1 Pure collinear contribution 

The collinear contribution is simple to assess since one can, to the NLL accuracy we seek, 
consider collinear radiation as being included in the evolution of the incoming and outgoing 
hard quark jets. In order to derive the NLL pure collinear contribution it proves useful 
to first consider the observable as defined in Eq. (|l.l|). Since one is dealing with soft 
and/or collinear gluon emission alone, we are looking at a tiny deviation from the Born 
configuration. Thus the sum over current hemisphere emissions, on the RHS of Eq. (|l.lj), 
includes a contribution from the transverse momentum of the outgoing current quark. To 
work in terms of secondary emissions alone, one uses conservation of transverse momentum 
to write Eq. ( p!j] ) as: 

Qt = -J2 > (2-1) 

where the sum now runs over all final-state particles in the remnant hemisphere. In the 
collinear region there is an important simplification, in that all these emissions can, to our 
accuracy, be attributed to the showering of the incoming quark. Note that since one is now 
inclusive over current-hemisphere emissions, we can neglect the collinear evolution of the 
outgoing quark. This will correct our resummed result by a factor of relative order as, but 
not enter into the NLL form-factor we aim to compute. 

The next step is to consider multiple collinear gluon branchings on the incoming hard 
leg. In this region the squared matrix-element can be approximated to NLL accuracy by 
a product in A'^ space of individual gluon emissions from the hard incoming quark, where 
N is the moment variable conjugate to Bjorken x. Taking first just soft and collinear 
emissions^, where we can just as well work in x space, one can write: 



1 da 



<To dQ 



where dPn is the differential n gluon emission probability and we introduced the vector 
Pt, which is the vectorial sum of transverse momenta of particles in TCji. Note that for 
purely soft and collinear emissions the parton distribution function (pdf) cancels with that 
in (To, the Born cross-section. This is not correct in the hard-collinear region and we will 
re-introduce the pdf while considering those emissions. For the final result, we shall also 



*The subsequent extension to hard and collinear emissions will be straightforward. 
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normalise to the cross-section including up to 0{a1) corrections, rather than merely the 
Born cross-section. 

We first compute the integrated quantity: 

-^{Q,Qt) = - r -^dQ'^^Y. I dPnO{Qt-pt)d'pt6' \pt+ • (2-3) 



One can then express: 

[pt + E = 1 n ^''^^ ■ (2-4) 

Having achieved our aim of factorising the delta function constraint into a product of 
individual gluon contributions, we integrate over pt which reduces the above to: 



R 



where in writing the above we made use of J^^ d6 ex.p (i b pt cos 9) = 27r Jo{bpt), and 
uJi{u) = u' Jo{u')du' with Jq and Ji being the zeroth and first order Bessel functions. 

The emission probability dPn factorises for soft and coUinear emissions into an essen- 
tially classical independent-emission pattern: 

dPn = -7 Cf ■ To^ar/i , (2.6) 

where kt^i and rji refer to the transverse momenta and rapidity of the i^^ emission with 
respect to the incoming quark direction . Since the incoming quark is anti-parallel to the 
photon axis, in the Breit frame, the kt immediately above is the same to NLL accuracy 
as that measured with respect to the photon axis and we do not distinguish the two. The 
coupling as is defined in the CMW scheme We note that the rapidity integration is 
bounded by = at large angles to the incoming quark since we are considering emissions 
in the remnant hemisphere alone. 

Summing over all emissions in the remnant hemisphere using the factorised forms ( |2.^ ) 
and ( ^.6| ) and inserting virtual corrections according to the emission pattern ( |2.6D (with 
an additional factor (—1)") we arrive at the resummed soft and collinear contribution to 
a{Q,Qt)/ao: 

J QtJiibQt) eM-Rsc{b)]db, (2.7) 

where Rscib) is the "radiator" accounting for soft and collinear emissions by the incoming 
quark: 

i?3,(6) = I ^ dvasik^) (e^^-^^' - l) . (2.8) 
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We remind the reader that one needs to correct the above expression to obtain single- 
logarithms arising from hard collinear radiation as well as those that arise in the large-angle 
region from soft emissions, which we shall do presently. 

Let us for the moment concentrate on the quantity Rsc{b), which is the most important 
piece of the result since it contains the leading (double) logarithms. Integrating over the 
polar angle variable in Eq. ( |2.8[ ) we obtain: 

2Cf f dkf 2\ 



Rscib) = / -j^driasikf) (Mbkt) - 1) . (2.9) 

TT J kt 

Further to next-to-leading or single logarithmic accuracy, it suffices to make the substitu- 
tion {jQ{bkt) — 1) ^ —9 {kt — 2e~^^ /b) and arrive at: 

i?sc(6) = ^ r a,{k^)^ln9., -b = be^-/2, (2.10) 
vr Ji/b kt kt 

where we performed the rapidity integration. 

Next we extend the soft-collinear result above to the full collinear one by including 
hard emissions. As is easy to show (see e.g. Q), hard emissions collinear to the incoming 
quark lead to a modification of the factorisation scale in the pdfs to the scale 1/6^, 
q{x, fi^) ^ q{x,l/lP). A remnant of this change of scale is the replacement of Q in Eq. ( |2.1C| ) 
by Qe~^/^. Thus the extension of the soft-collinear result, Eq. (|2.7| ), to the full collinear 
one is: 

^ ^' q{x,l/b^)QtJiibQt) exp[-R(bQ)]db, (2.11) 



g(x,/i2) 

with: 



2Cf r n2^dkt A„ Q 3 



vr 



1/6 



R{hQ) = ^ / In^-- . (2.12 



2.2 Non-global corrections 



We also have to include the effects of soft emissions at large angles. Thus far we have identi- 
fied remnant emissions as those belonging to the incoming quark while current hemisphere 
emissions (over which we claimed to be inclusive) are associated to the struck final-state 
quark. As is the case for single-hemisphere observables, this identification is not correct at 
single-log level due to correlations between soft emissions at large angles [11, 12|. Thus the 



remnant hemisphere distribution is affected at SL level by soft gluons at large angles to the 
current quark, but still in the current hemisphere, emitting into the remnant hemisphere. 

The computation of this piece has been carried out, in the large Nc approximation^ 
and is universal for all observables having a linear dependence on kt of soft large-angle 
emissions. We label this "non-global" piece S{bQ), which can be parameterised as |11|: 

SibQ)^ exp {-C^C^^ (iTM^) '1 ' ^'-''^ 



^Strictly speaking, the full result has been computed at 0{ag) and the large Nc approximation starts 
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where: 

*(^^) = -j^l^{'i--'^o^s{Q'')l3oHbQ)) , (2.14) 

with A = 0.85 Ca, B = 0.86 Ca and C = 1.33. 

Thus our final form for the resummed result to NLL accuracy can be expressed as: 

— a(Q,Qt) = — i-^r f q{xMb')S{bQ)e-'^^^'^^QtJi{bQt)db, (2.15) 

where we chose ji = Q. The above result now incorporates all the sources of logarithmic 
enhancements to NLL or single-log accuracy, specifically soft and hard collinear emissions 
and soft emissions at large angles. The result for the NLL radiator R(bQ) is explicitly 
given in appendix ^. In the following section, we shall take the b space result above and 
convert it to one valid in Qt space, over the range of Qt values that interest us. 

3. Result in Qt space 

We start by noting that one commonly used method to derive a Qt space result from 
the b space form is simply to evaluate the complete b integral in Eq. ( p. 15 ) numerically. 



This method is not without several well-documented shortcomings |l^, that are usually 
circumvented by "reasonable" prescriptions that are not derived from first principles of 
QCD. 

For instance the b integral stretches from to cxd but the function R{b Q) has a Landau 
pole singularity at bQ = exp{l/ (2/?o Oig)}, which means it is perturbatively undefined for 
large b values. To get around this problem one introduces a parameter b* and substitutes 



6 — > 6* = Vyl + ^^/^Um 111- This ensures one never evaluates R(bQ) or the structure 
functions at scales larger than some cut-off bum, whose value is adjustable. Additionally, 
we smear the b space result with a non-perturbative Gaussian function that is also not 



obtained from first principles, but typically through fits to data sets |20, 21, 
These prescriptions are needed in order to do the b integral and obtain a result for finite Qt, 
even at relatively large Qt values where one might expect to trust perturbative predictions 
and where additional non-perturbative parameters or ad-hoc inputs should not play any 
significant role. Moreover the matching to fixed-order Qt space results is complicated by 
not having an analytical resummed result in Qt space. 

In what follows we provide an analytical Qt space result which is valid for use over a 
large range of Qt values and represents a clean extraction of the next-to-leading logarithmic 
resummed Qt space result, from the b integral. The price we pay for not evaluating the 
complete b integral in detail, is a formal divergence at small Qt, which one can anticipate 
quite generally through considerations based on the work of Parisi and Petronzio ||l^ . Thus 
we cannot use our result at very small Qt values since in this region our answer is no longer 
a good approximation to the b integral. For quantitative studies however, with the specified 
Qt range over which data is available, our formula is valid for use as it stands. The region 
over which we start to see a problem with our approximation occurs at Qt values that are 
too small to study accurately via perturbative methods and in any case below the lowest 
Qt data. 
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To obtain a resummed result in Qt space we expand the function R{bQ) in Eq. ( 2.15| ), 
about the point 6(= hQt) = 2e~^^ to obtain: 

R{bQ) = RiQ/Qt) + R'iQ/Qt) (-fE-ln2 + \nb)+---, (3.1) 



where we have used R'{bQ) = dR{bQ)/dlnb and neglected R" and higher derivatives as 
they contribute only to subleading (below single-log) accuracy. The non-global function 
S and the b dependent parton distributions are straightforwardly expressed in Qt space 
with the substitution b — > 1/Qt, since logarithmic derivatives of these single-log functions, 
analogous to R' above, only contribute at subleading accuracy. 

Using the Taylor expansion for R above, one can cast the result ( p. 15 ) as: 



La{Q,Qt) = i(^1^5(Q/Qi)e-«(«/'3')+('"2-^-)«'(«/«') T d6Ji(6) ^^'W/'?') . (3.2) 
c^o q{x, Q^) Jo 

Integrating over the dimensionless quantity b and incorporating, as a factor, the 0{as) 
constant pieces (see appendix ^) we can write the result as: 

X SiQ/Qt) e-^(Q/Q0-7.«'(Q/Q^) ^|j--|^;/^| , (3.3) 

where Cq and Ci are matrices in flavour space of coefficient functions (see appendix 
0). While the Co terms are merely proportional to delta functions the Ci pieces are 
important to correct the soft-collinear resummed result for hard real and virtual emissions 
at the leading 0{as) accuracy. They are straightforward to compute and are presented in 
appendix 

We immediately note that the above result diverges at R' = 2, an entirely expected 
feature. The reason for this divergence (encountered also in the Drell-Yan Qt resummations 
and the jet broadening in DIS ||^]) is merely the fact that at very small Qt the result one 
obtains is not described by exponentiation of the leading-order result, which is essentially 
the form we have derived above |l^. As — > the mechanism of vectorial cancellation 
between emissions of formally arbitrary hardness, takes over from the Sudakov suppression 
of soft and collinear radiation, as the dominant mechanism for producing a small Qt- 
However the divergence does not play a major role for phenomenological purposes, since 
over the values of Qt we intend to study we are sufficiently away from the point R' = 2. 
This will be further elaborated in the subsection below. 

3.1 Position and impact of the divergence 

As we mentioned above, for the particular case at hand, the divergence occurs at rather 
small values of Qt for the Q values^ of interest to us. The corresponding Qt values, at and 
near the divergence, fall in a region that is either neglected for phenomenological purposes 



^There are data in the range 17GeV < Q < 116 GeV. 
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or modeled with the introduction of non-perturbative parameters, since one expects non- 
perturbative effects to be large here. In the Qt region where the divergence does not have 
any significant numerical impact, we still probe small enough Qt to test the perturbative 
resummation and non-perturbative corrections, as is our aim. 

To be precise, the divergence occurs at R' = 2. In terms of the variable Qt, using the 
expression for R' given in appendix Q this results in: 

Qt = Q exp ( ^ \ , (3.4) 

which for the illustrative value of Q = 90 GeV gives Qt = 0.52 GeV, with as = 0.118. 
Since this region is in any case perturbatively unsafe, being quite close to the QCD scale, 
we do not expect to obtain sensible results with the perturbative methods we have used. 
However we can safely study Qt values of a few GeV without worrying about the impact 
of the formal divergence at R' = 2. 

We can quantify this statement as follows: in the region where R' = 2, there is a 
complete breakdown of the hierarchy between leading, next-to-leading, etc. logarithms. In 
order to determine up to which value of R' one can use the usual hierarchy, where N"LL 
terms are suppressed by a" with respect to LL terms, one can follow the procedure outlined 
in 0. 

From those arguments one can infer that terms that are formally NNLL contribute a 
correction that is of the same order as the terms that one keeps in the NLL resummed 
result, in the region where 2 — R' = The NNLL terms contribute at relative 0{as) 

when i?' = 1 or more. Thus for i?' < 1 we can safely use our resummed Qt space formula 
since omitted NNLL terms contribute as usual, at relative 0{as)- 

We can work out the position of both points in terms of Qt for a given Q. For Q = 90 
GeV we obtain that the critical value, where all terms in the formal hierarchy are in fact 
of the same order, is Qt = 0.68 GeV and that the region where the usual hierarchy is 
respected, is Qt > 1.5 GeV. This still allows the full range of available data to be safely 
probed, including the lowest measured Qt bins. 



4. Matching to fixed-order 

Having obtained the NLL perturbative estimate we now need to combine it with the exact 
C(a^) perturbative result to obtain accurate predictions over the entire range where data 
exist. We follow the matching prescription known as M2 matching introduced in Ref. Q. 
Here, the final result is given by: 

ar + as (41) - 41)) + al (af) _ ^2)^ ^^q^ ^ (4^^) 

where Os = a<j/(27r) and ai^^'^"^^ denote the coefficients of the resummed result, ar, ex- 
panded to first and second order in ag respectively, while al^^'^"^^ are the corresponding 



coefficients obtained from fixed-order Monte Carlo programs such as DISENT |14]. The 
above matching formula adds the resummed and exact results and subtracts the double- 
counted terms (up to 0{a'^)) that are included in the resummation. 



-9- 




Figure 1: Comparison between DISENT IQ, matched resummed and pure resummed differential 
distributions. MRST NLO pdfs have been used with = 0.1205 [||. 

Note that terms such as a'^ln{Q/Qt) that are formally subleading and hence not in- 
cluded in the resummation, are present in the piece (^a^^ — cr^^^^ of Eq. ([4.1|). This 
piece is divergent as — > and thus we multiply it by the resummed form-factor 
T,{Q,Qt) = q{x,Qf)/q{x,Q^)S{Q/Qt)ex.p{—R{Q/Qt)}, to ensure sensible behaviour at 
small Qt- This procedure is obviously ad-hoc but only affects the result at subleading 
logarithmic accuracy, which is in any case beyond our quantitative control. 

Another point that needs to be re-emphasised is that the factor T, as we use it here, is 
just an approximation to the resummed result given by a full evaluation of the b integral. 
The approximation is intended for use (and is valid to NLL accuracy) only sufficiently away 
from R' = 2, the position of the divergence. As we explained in the [last| section, this covers 
the range over which data exist and over which we can make reasonable comparisons. 

5. Results 

The aim of this section is to display the results obtained for a next-to-leading logarithmic 



resummed prediction matched to next-to-leading order predictions from DISENT |14|. 
Additionally we comment on the role that might be played by a non-perturbative Gaussian 
smearing function and that small-x effects may effectively give an enhanced smearing of 
the spectrum leading to a broader prediction than the one provided here as was observed 
also in the SIDIS case |, |]. 

In Fig. |l], we display the fixed order NLO results along with the pure resummed and 
matched results for two values of the hard scale Q and corresponding Bjorken x values. The 



choice of the x and Q values correspond to bins where preliminary data already exist [13|. 
There one notes the divergent NLO result and the correct small Qt behaviour as given by 
the resummation as well as the role of matching in the high Qt tail of the distribution. We 
also point out that the role of the non-global term is limited to a few percent effect after the 
matching to fixed order has been performed. Thus missing uncalculated single-logarithmic 
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Figure 2: Non-perturbative smearing of the resummed prediction with a Gaussian function in kt, 
exp[— g/c^]. Two different choices of the smearing parameter g are illustrated. 



terms in the non-global piece S that are suppressed as are not expected to change 

our quantitative conclusions. 

We also present here the effect of smearing or convoluting our resummation with a 
b space or equivalent kt space Gaussian function representing non-perturbative effects. 
Such intrinsic kt effects have been the subject of much phenomenological investigation in 
Drell-Yan like processes |5|, |2^, ^ 27 1. The non-perturbative smearing function we choose 
has the simple form F^^{kt) = exp[—gkt] and we perform a two-dimensional convolution, 
/ (Pkt g/irF^^ (\Qt — kt\^ I^^t) with the resummed distribution I{kt) = l/ada/dk^, to 
obtain the smeared result as a function of Qt- We carry out the smearing with different 
values of g, two of which are illustrated in Fig. ^. The value ol g = 0.5 GeV~^ corresponds 
to a reasonable mean intrinsic kt value of 1.25 GeV. We also illustrate the effect of using 
a smaller g value of 0.1 GeV~^, which leads as shown to a broader Qt distribution. Ideally 
one would need to compare our predictions with the data at different x values in order to 
ascertain whether one sees any visible broadening of the NLL resummed spectra at smaller 
X, such as that mimicked by a change in the smearing parameter g. If this is the case one 
may investigate the dependence of the smearing parameter g on Bjorken x, a step we leave 
to forthcoming work p5| ]. 

We should point out that in convoluting the resummed prediction with the non- 
perturbative Gaussian kt distribution it was necessary to provide an extrapolation of our 
resummed results for I{kt) down to kt = 0. We have chosen this extrapolation as in 
Ref. 1 18 1 so as not to modify the NLL resummed result as well as to obtain the correct 
limiting behaviour of the integrated cross-section oc /cj , as determined by evaluating the h 
integral for the integrated quantity (j{kt), in the limit kt 0. 

Thus we substitute for kt an effective variable (fc^ )^ = kt + Qq exp {—kt/Qo) and make 
the replacement a{kt) — > f (A;^) (l — exp (— a/c|)), where Qo and a should be chosen so as 
not to modify the resummed result in the range where it is valid, as in Ref. ll^. These 
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additional parameters should only modify the very low (non-perturbative) Qt region where 
our NLL result is in any case not valid as it stands. One can think of the parameters a and 
Qq as being non-perturbative inputs that can be varied alongside the smearing parameter 
g to obtain a good fit to data in the lowest Qt region. 

For the plot in Fig. |2|we chose to take Qq = 1.2 GeV and a = I/Qq as these choices do 
not impact our resummed results over most of their range of validity. The dominant impact 
in the very low Qt region, beyond the control of NLL resummation, is in fact that of the 
smearing function exp (—gkt), for which different choices have been already mentioned. 



6. Conclusions 

We have introduced here a DIS variable that, as we have explained, has a very simple re- 
lationship to vector-boson and Higgs Qt spectra at hadron colliders. The aim of doing this 
has been to use HERA data to compare resummed theoretical predictions with experiment, 
keeping an eye on comparions at lower x values. This complements the extensive studies of 
DIS event shapes that have been carried out thus far which employed the standard resum- 
mation formalism (and non-global logarithms), supplemented by l/Q power corrections 



| 10f| . We recall that the program of comparing DIS event shapes to the data was quite 
successful without any special role visible for small-x effects . 

Given however that only moderately small x values, x ~ 10~^, can be reasonably 
accessed in these studies, it is clearly better to choose a variable that is potentially more 
sensitive to small- re dynamics than event-shape variables, in order to determine the role of 
these effects. We expect such a variable to be the Qt spectrum we have defined here, where 
it will be interesting to establish if a small-x enhanced broadening of the resummation we 
presented, is indeed visible in the data. This was apparently the case in SIDIS studies Q 
and, if present, we expect these effects to manifest themselves for our observable too. Given 
the simple relationship of our results to those for Drell-Yan like observables it should then 
be easy to extrapolate our conclusions to hadron collider studies, where it is important to 
reach a firm conclusion on the issue of small- a; broadening. 

At present we have preliminary data for our observable from HERA and we await the 
data in its final form before making detailed comparisons and drawing phenomenological 
conclusions on this issue. This will be the subject of forthcoming work. 
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A. Leading order result 



We report briefiy below the leading-order result for crr{Qt,Q), where r denotes the re- 
summed result, neglecting terms that vanish as Qt — > 0. In line with the procedure of 



Ref . M , the leading-order cross-section for events with 



< Qt normalised to the 
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Born cross-section is^: 



+ 



—2 \ -J- I Q 



'-2 \ -7- I ^ Q 



. (A.l) 



where y is the usual DIS Bjorken variable, Tf = TrY^^-c^ and in the MS scheme we have: 



21 



- (12^' - loe' + 1) 



1+^ 

g(2e-l) 
2(1-0 



1 

In— + 



, (A.2) 



+ 



■^.,2 (C, ^) = -2[e^ + (1 - 0'] In ^ + 2^(1 - 0+ 



+ {{2i- 1){-Qe + 6^ - 1) - [e' + (1 - 0'] In ^ 



^(2^-1), (A.3) 



J^,,i = 2^(2^-1) 0(2^-1), 



(A.4) 



= 4^(1 - 0(2? - 1)^ (2? - 1) . (A.5) 
The (2n/ + 1) x 1 matrices Cq and Ci are defined such that their transposes are given 



by: 



and: 



elS{l-z) 



V 



Cj{z) 



Cf el {j'qAz, 1) - ^^((_yyl J^qx{z, 1)} ^ 
Cf el 1^9,2(2, 1) - i+(f^-^g,L(2, 1) 

•^.,l(^,1)} J 



(A.6) 



(A.7) 



'^We compute both the F2 and Fl contributions indicated by the subscripts 2 and L 
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B. The radiator 



The radiator is given by: 



R{Q/Qt) = LgiiasL) + g2{asL) , 



(B.l) 



where L = ln{Q/Qt) and in the MS scheme: 



Cf_ 

7r/3o 



1 + 



ln(l-2A) 
2A 



(B.2) 



47r/?o 



with A = as{Q'^)(3oL and: 



ln(l - 2A) + 



2A 



1 - 2A 



27r/33 



2A + ln(l-2A) 1 

1-2A 2 ^ ^ 



, (B.3) 



K = Ca 



67_^ 
18 6 



-n 



lie A -2nf 
12^ 



/3i 



17Cl-5CAnf -SCpUf 
24^ 



(B.4) 

In obtaining this, we used Eq.( p.l2 ) and the 2-loop QCD (3 function to replace the scale of 
as with and moved from the CMW |17| scheme to the MS scheme (see e.g. Ref. ||7|). 
The derivative of the radiator with respect to ln6 at = Q/Qt is given by: 

2Cf a 



R'{Q/Qt 



7r/3o 1 - 2A 



(B.5) 



The expansion of the resummed result (Eq.(|3.3D) to 0{as) and 0{a'^), which is needed 
in Eq. (O) , yields: 



(1)/ n r^r C^®P^q{x^^ Ci®q(x,Q^) 



(B.6) 



= G22L' + + ^Gl^L^ - 2 ((Gil + 27r/3o)L2 + G12L' 



g(x,Q2) 



^2\ 



+ g(x,Q2) (GiiL + Gi2L)+2 L, (B.7) 

where P is the matrix of leading order splitting functions and the coefficients Gmn are 
given in table |l[ One can clearly see that the expansion of the resummed result to 0{as) 
reproduces the leading order result given by Eq. (A.l). 
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Gl2 


-2Cf 


Gil 


3Cf 


G23 




G22 


-^CfCa + I C2 - 2 CfK + 6 tt CfPo 



Table 1: The coefficient Gnm that enter the fixed-order expansion of the resummed resuh. 
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